A family of density expansions for Levy-type processes with default 
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Abstract 

We consider a defaultable asset whose risk-neutral pricing dynamics are described by an exponential 
Levy-type martingale subject to default. This class of models allows for local volatility, local default 
intensity, and a locally dependent Levy meas u re. G eneralizing and extending the novel adjoint expansion 
technique of iPagliarani. Pascucci. and Rieal (|2013l ). we derive a family of asymptotic expansions for 
the transition density of the underlying as well as for European-style option prices and defaultable 
bond prices. For the density expansion, we also provide error bounds for the truncated asymptotic 
series. Additionally, for pure diffusion processes, we derive an asymptotic expansion for the implied 
volatility induced by European calls/puts. Our method is numerically efficient; approximate transition 
densities and European option prices are comp uted via Fourier transfo rms; approximate bond prices are 
computed as finite series. Additionally, as in 



Pagliarani et al. 



( 20131 ). for models with Gaussian- type 



jumps, approximate option prices can be computed in closed form. Numerical examples confirming the 
effectiveness and versatility of our method are provided, as is sample Mathematica code. 

Keywords: Local volatility; Levy-type process; Asymptotic expansion; Pseudo-differential calculus; De- 
faultable asset 



1 Introduction and literature review 

A local volatility model is a model in which the volatility at of an asset X is a function of time t and the 
present level oi X. That is, at = a{t, Xt). Amo ng lo c al vol atility models, perhaps the most vifcU-known is the 
constant elasticity of variance (CEV) model of lCoxl (jl975h . One advantage of local volatility models is that 
transition densities of the underlying - a,s well as Europ ean option prices - are often available in closed-form 
as infinite series of special functions (see lLinetskvl (J2007I ) and references therein). Another advantage of local 
volatility models is that, for models whose transition density is not available in closed f orm, a ccurate density 
and opt ion price app roximations are readily available (see, IPagliarani and Pascuccil (|2011l ). for example). 
Finally, iDupird (jl994[ ) shows that one can always find a local volatility function cr(i, x) that fits the market's 
implied volatility surface exactly. Thus, local volatility models are quite flexible. 
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Despite the above advantages, local volatility models do suffer some shortcomings. Most notably, local 
volatility models do not allow fo r the underlying to experience jumps, the need for which is well-documented 
in literature (see lErakeiJ (|2004r ) and references therei n). Recently, there has been mu ch interest in combin- 



ing local volatility models and models with ju mps 



Andersen and Andreasen 



( 20001 ). for example, discuss 



extensions of the implied diffusio n approach oflDupird (119941) t o asse t processes with Poisson jumps (i.e., 
jumps with finite activity). And iBenhamou. Gobet. and Miril (|2009l ) derive analytically tractable option 
pricing approximations for models that include local volatility and a Poisson jump process. Their approach 
relies on asymptotic expansions aroun d small diffusion and small jump frequency/size limits. More recently, 
Pagharani. Pascucci. and Rigal (120131) consider genera l local v olatility models with in dependent Levy jumps 



Benhamou et al 



(l2009l) . 



Pagliarani et al 



(|2013f ) make no small jump 



(possibly infinite activity). Unlike, 
intensity/size assumption. Rather the authors construct an approximated solution by expanding the local 
volatility function as a power series. While all of the methods described in this paragraph allow for local 
volatility and independent jumps, none of these methods allow for state- dependent jumps. 



Stochastic jump-intensity was recently identified as an important feature of equity models (see lChristoffersen. Jacobs, and ( 



(|2009l )). A locally dependent Levy measure allows for this possibility. Recently, two different approaches have 



Mendoza-Arriaga. Carr. and Linetskv 



been t aken to modeling assets with locally-dependent jump measures 
(|2010l ) time-change a local volatility model with a Levy subordinator. In addition to admitting exact option- 
pricing formulas, the subordination technique results in a locally-dependent Levy measure. iLorig] (J2012r ) 
considers another class of models that allow for state-dependent jumps. The author builds a Levy-type pro- 
cesses with local volatility, local default intensity, and a local Levy measure by considering state-dependent 
perturbations around a constant coefficient Levy process. In addition to pricing formula, the author provides 
an exact expansion for the induced implied volatility surface. 

In this paper, we consider scalar Levy- type processes with regular coef ficients, which naturally includes all 



Pagliarani et al 



(|2013l ). we derive 



the models mentioned above. Generalizing and extending the methods of|] 

a family of asymptotic expansions for the transition densities of these processes, as well as for European-style 

derivative prices and defaultable bond prices. The key contributions of this manuscript are as follows: 



Pagliarani et al. 



Weallow for a locally-dependent Levy measure and local default intensity, whereas 
(|2013[ ) consider a locally independent Levy measure and do not allow for the possibility of default. 
A state-dependent Levy measure is an important feature because it allows for incorporating local 
depende nce into infinite activity Levy m odels that have no diffusion component, such as Varianc e 
Gamma (Madan. Carr. and Chang ( 19981) ) and CGMY/Kobol ( Bovarchenko and Levendorskiil (|2000l ); 



Carr. Geman. Madan. and Yoij (12002)) 



• Unlike 



Benhamou et al. 



ij 



(|2009l ). we make no small diffusion or small jump size/intensity assumption. 



Our formulae are valid for any Levy type process with smooth and bounded coefficients, independent 
of the relative size of the coefficients. 
Whereas ' 



iPagliarani et al.l (|2013l ) expand the local volatility and drift functions as a Taylor series about 
an arbitrary point, i.e. f{x) = ^„ a„(a;— a;)", in order to achieve their approximation result, we expand 
the local volatility, drift, killing rate and Levy measure in an arbitrary basis, i.e. f{x) — J2n (^nBn(x). 
This is advantageous because the Taylor series typically converges only locally, whereas other choices 



of the basis functions _B„ may provide global convergence in suitable functional spaces. 
Using techniques from pseudo-differential calculus, we provide explicit formulae for the Fourier trans- 
form of every term in the transition density and option-pricing expansions. In the case of state de- 
pendent Gaussian jumps the respective inverse Fourier transforms can be explicitly computed, thus 
providing closed form approximations for densities and prices. In the general case, the density and 
pricing approximations can be computed quickly and easily as inverse Fourier transforms. Additionally, 
when considering defaultable bonds, approximate prices are computed as a finite sum; no numerical 
integration is required even in the general case. 
For models with Gaussian-type jumps, we provide poin t wise error estimates for transition densities. 



Thus, we extend the previous results of iPagliarani et al.l (|2013l ) where only the purely diffusive case is 



considered. Additionally, our error estimates allow for jumps with locally dependent mean, variance 
and intensity. Thus, f or models with Gaussian-type jumps, our results also extend the results of 



Benhamou et al. 



(|2009f ). where only the case of a constant Levy measure is considered. 
• For local volatility models with no jumps and no possibility of default, we provide an asymptotic 
expansion for implied volatilities corresponding to European calls/puts. As with approximate bond 
prices, approximate implied volatilities are computed as a finite sum of simple functions; no numerical 
integration is required. 

The rest of this paper proceeds as follows. In Section, ^ we introduce a general class of exponential 
Levy-type models with locally-dependent volatility, default intensity and Levy measure. We also describe 
our modeling assumptions. Next, in Section|3l we introduce the European option-pricing problem and derive 
a partial integro-differential equation (FIDE) for the price of an option. In Section |4] we derive a formal 
asymptotic expansion (in fact, a family of asymptotic expansions) for the function that solves the option 
pricing FIDE (Theorem |8]). Next, in Section [5l we provide rigorous error estimates for our asymptotic 
expansion for models with Gaussian- type jumps (Theorem 1 14|). In Section [6] we derive an implied volatility 
expansion for pure diffusion models (i.e., models with no default and no jumps). This expansion involves no 
special functions and no integrals and can therefore be computed extremely quickly. Lastly, in Section [71 we 
provide numerical examples that illustrate the effectiveness and versatility of our methods. Technical proofs 
are provided in the Appendix. Some concluding remarks are given in Section [8l 

2 General Levy-type exponential martingales 

For simplicity, we assume a frictionless market, no arbitrage, zero interest rates and no dividends. Our 
results can easily be extended to include locally dependent interest rates and dividends. We take, as 
given, an equivalent martingale measure Q, chosen by the market on a complete filtered probability space 
(ri, 3^, {3^t, t > 0}, Q) satisfying the usual hypothesis of completeness and right continuity. The filtration J't 
represents the history of the market. All stochastic processes defined below live on this probability space and 
all expectations are taken with respect to Q. We consider a defaultable asset S whose risk-neutral dynamics 



are given by 



St 



kc>t}e''\ 



dXt== fi{t,Xt)dt + cr{t,Xt)dWt+ / dNt{t,Xt^,dz)z, 
dNt(t,Xt^,dz) ^dNt(t,Xt-.,dz) - v{t,Xt^,dz)dt, 



C = inf U > 



7(s,X,)ds>£ 



(1) 



Here, X is a Levy-type process with local drift function /j,(t, x), local volatility function (j{t, a;) > and state- 
dependent Levy measure v{t, x, dz). We shall denote by 3^^ the filtration generated by X. The random vari- 
able £ ~ Exp(l) has an exponential distribution and is independent of X. Note tha t C. which represents the 
defaul t time of S, is constructed here trough the so-called canonical construction (see 



Bielecki and Rutkowski 



pOOll )). and is the first arrival time of a doubly stochastic Poisson process with local intensity function 



^(t,x) > 0. This way of modeling default is also considered in a local volatility setting in 



Carr and Linetskv 



(J2006l ): iLinetskvl (|2006f ). and for exponential Levy models in ICapponi et al.l (|2013r ). Since C is not Jj 



F^ 



measurable we introduce the filtration 3"^ = o" ({C ^ s) , s < i} in order to keep track of the event {C, < t}. 
The filtration of a market observer, then, is J't = 'J^ V 3^f . 

In the absence of arbitrage, S must be an 3^t-martingale. Thus, the drift fi{t,x) is fixed by a{t,x), 
v{t, X, dz) and 7(i, x) in order to satisfy the martingale condition |^ 

IJ,{t,x) — j{t,x) — a{t,x) — / i'{t,x,dz){e'^ — 1 — z), a{t,x) := -a'^{t,x). (2) 

In order for equation ([Ij and subsequent computations to make sense, we will assume the following: 

1. Smoothness. The coefficients are continuous in t and smooth in x: a(-, •), 7(-, •), lyO, •, dz) G C^'^^fM) 

i 

2. Boundedness. The coefficients a and 7 are bounded. There exists a Levy measure 

v{dz) := sup i/{t,x,dz) 

{t,x)eR+xR 

such that 

z>((iz) min(l, z^) < 00, / z>((iz)e^ < 00, / P{dz)\z\< 00. (3) 

J\z\>l "'|z|>l 

Note that these conditions imply the existence of a positive constant p, < co such that |/i(i, x)\ < p, for 
all (i,x). 

3. Lipschitz Continuity. There exists a constant C such that 



\lJ.(t,xi) - fi(t,X2)\'^ + |cr(t,xi) -CT(i,X2)P + I z^ ^/v{t,xi,dz)- ^Jv{t,X2,dz) < C\ 
for any t G K+ and xi, 2:2 G M. 



Xi - X2 



^ We provide a derivation of the martingale condition in Section [3l Remark [2l below. 
^ This notation for 1/ means; v(-, -,-8) G C'''°°(]R) for any Borel set B. 



4. Existence of a density. Either a(-, •) > or for every a; e M there exists a{x) £ (0, 2) such that 

liminfe"(=")~2 / z^iy(x,dz) > 0. 

These conditions guarantee that the distri bution of X t has density with respect to the Lebesgue measure 
(in fact, a C°° density) for all t > 0. See lSatol (|l999l ) Proposition 28.3. 



Toget her, conditions [2] and 



3] ensu re that the SDE in ([T]) has a unique strong solution up to the default time 



C (see lOksendal and SulemI (|2005l ). Theorem 1.19). We include Assumption 01 because, in what follows, we 
speak of densities. Nevertheless, the existence of the density is not strictly necessary in our analysis. Indeed, 
since our formulae are carried out in Fourier space, all of our computations are still formally correct even 
when dealing with distributions that are not absolutely continuous with respect to the Lebesgue measure. 

We will relax some of these Assumptions for the numerical examples provided in Section [71 Even without 
the above Assumptions in force, our numerical results indicate that our approximation techniques work well. 



3 Option pricing 

We consider a European derivative expiring at time T with payoff H{St) and we denote by V its no-arbitrage 
price. For convenience, we introduce 

h{x):=H{e') and K := H{Q). 

Proposition 1. The price Vt is given by 



Vt=K + ^^yt}^ 



- S;^ lis,X,)ds 



{h{XT) - K) \Xt 



t<T. 



(4) 



The proof can be found in Section 2.2 of 



m 



Linetskvl (|2006f) . Because our notation differs from that of 



Linetskvl (|2Q06r ). and because the proof is short, for the reader's convenience, we provide a derivation of 
Proposition [H here. 

Proof. Using risk-neutral pricing, the value Vt of the derivative at time t is given by the conditional expec- 
tation of the option payoff 

Vt=W.[H(STWt] 

= E [h{XT)l{c>T}\'^t] + KE [l^^^T] \^t] 
= E [h{XT)l{c>T}\3't] +K-KE [I{^>T} \1t] 
^K + I{c>t}]E [e- i^ ^i^^^s)d^ {h{XT) - K) \'jf 
^K + l{c>t}E \e- ^" ^(-■^=)'^^ {h{XT) - K) \Xt 



where we have used Corollary 7.3.4.2 from lJeanblanc. Yor. and Chesnevi ()2009l) to write 



D 



Remark 2. By Proposition [T] with K ^ and h{x) = e^, we have that the martingale condition St = 
E[ST\3^t] is equivalent to 

Therefore, we see that S" is a martingale if and only if the process exp ( — /n 7 (s, Xg) ds + Xtj is a martingale. 

The drift condition ^ follows by applying Ito's formula to the process exp (— L ^ (s, Xs) ds + XA and 
setting the drift term to zero. 

From (J4|) one sees that, in order to compute the price of an option, we must evaluate functions of the 
formlfl 



v{t,x) :=E 



(5) 



g-/-^(s,x.)ds/^(^^)|^^^^^J ^ J h{y)p{t,x;T,y)dy 

Here, p{t, x; T, y) is a Feynman-Kac (FK) transition density, which represents the transition density of log 4j 

pit,x]T,y)dy = Q[\ogSTedy\logSt^x], x,y€R, < t < T < oo. 

Note that the FK transition density is not a probability density since (due to the possibility that St = 0) 
we have 

/ Pit,x]T,y)dy < 1. 

Jr 

It is well-known that the function v{t, a;), defined by ^, satisfies the Kolmogorov backward equation 

(dt + A)v = 0, v{T, x) = h{x). (6) 

Here, the operator A is the generator of the FK semigroup "Pi.T- That is. 



Af(x) = lim 



yt,Tf{x) - fix) 

T-t 



(if strong limit exists) 



(7) 



For f E Cq (K.) , the limit (O exists and the operator A is given by 



Af{x) - 7(i, x){d,f{x) - fix)) + ait, x)idlfix) - dJix)) 

(e- - 1 - z)vit, x, dz)dJix) + f ifix + z)- fix) - zdJix))vit, x, dz). (8) 

From now on we define dom(yi) to be all functions / such that, with A given by (j8]), the derivatives and 
integrals in Afix) exist and are finite for all x e R. Note also that, if one sets h = 5y, then vit,x) becomes 
the FK density pit, x; T, y) since 

6yiz)pit, x; T, z)dz = pit, x; T, y). 

From the analytical point of view, pit, x; T, y) represents the fundamental solution of dt + A. 



^Notc: we can accommodate stochastic interest rates and dividends of the form rt = r{t,Xt) and qt = q(t,Xt) by simply 
making the change: ■y(t, x) — > 7(i, x) + r{t, x) and fi{t, x) — > fj,{t, Xt) + r(t, Xt) — q{t, Xt). 
^Here with logS we denote the process Xtl^^yt) ~ °oI{c<t} 



Remark 3. If S is the generator of a scalar Markov process and dom(S) contains S(]R), the Schwartz space 
of rapidly decaying functions on K.. Then S must have the following form: 



S/(x) = -l{x)fix) + ^i{x)^.J{x) + a{x)dif{x) + / i^{x, dz){f{x + z) - f{x) - I{|.|<i,}z9,/(a;)), (9) 



where 7 > 0, a > 0, z^ is a Levy measure for every x and R G [0, 00] (see iHohl (|l998f ). Proposition 2.10). If 
one enforces on S the drift and integrability conditions 1^ and ([3]), which are needed to ensure that S is 
a martingale, and allow setting R = 00, then the operators ([5]) and ^ coincide (in the time-homogeneous 
case). Thus, the class of models we consider in this paper encompasses all non- negative scalar Markov 
martingales that satisfy the regularity and boundedness conditions of Section [51 

Remark 4. In what follows we shall systematically make use of the language of pseudo-differential calculus. 
More precisely, let us denote by 

M^)=MO = ^e'^^, x,^eR, (10) 



the so-called oscillating exponential function. Then A can be characterized by its action on oscillating 
exponential functions. Indeed, we have 

Atp^{x) = (l){t,x,^)'il;^(x), 

where 

(bit, X, = 7(t, xm - 1) + a{t, x){^e - *?) (11) 

- / v{t, X, dz){e'' - 1 - z)ii + / v{t, X, dz){e'^'' - 1 - i^z), 

Jr Jr 

is called the symbol of A. Noting that 

e'^^==u(a;) = ^ —d"u{x) = u{x + z), 

for any analytic function u{x), we have 

{u{x + z) — u{x) — zdxu{x)) v{t, X, dz) = / i'{t, x, dz) ("e^ — 1 — zdx) u{x). 

Jr 

Then A can be represented as 

A^(j,{t,x,'D), V^-idx, 

since by ([TT]) 

(j){t, x, D) = 7(t, x)idx - 1) + a{t, x){dl - dx) 

v{t, X, dz){e^ - 1 - z)dx + / v{t, x, dz) (e^^=" - 1 - zd^) ■ 



Moreover, if coefScients a{t) , ^(t) , iy(t, dz) are independent of x, then we have the usual characterization of 
A as a, multiphcation by operator in the Fourier space: 

where 3^ and If~^ denote the (direct) Fourier and inverse Fourier transform operators respectively: 
3^1(0 - m ■■=^ f e-'^^fix)dx, 7-^f{x) = ^ / e^'^fiOdi. 



4 Density and option price expansions (a formal description) 

Our goal is to construct an approximate solution of Cauchy problem ([S]), which we assume has a unique 
classical solution. Sufficient conditions for the existence an d uniqueness of solutions of second order elliptic 



integro-differential equations are given in Theorem II. 3.1 of iGarroni and Menaldil (|l992l ). To find a solution 



of ([6]) , we assume that the symbol of A admits an expansion of the form 

oo 
(k{t,X,i) = Y,Br,{x)<|,r^{t,(,), (12) 

n=0 

where (f>nit,0 i^ of the form 

i^n{t,dz){e^ -l~z)i^+ I Vn{t,dz){e'''^ -l-izCl- (13) 



and {-Bn}n>o is some expansion basis with _B„ being an analytic function for each n > 0, and Bq — \. Note 
that (j)n{t,S,) is the symbol of an operator 

A„:=</)„(i,2)), D^-td^, (14) 

so that 

AnM^)^^n{t,OMx)- (15) 

Thus, formally the generator A can be written as follows 

oo 

■A = J2 Bn{x)An. (16) 

n=0 

Note also that Bn{x) is the symbol of the differential operator Bn{^id^) since 

Bn{-id^)MO = S„(a:)^,(e), n > 0. (17) 

Typically, An is the generator of a time-dependent Levy-type process (i.e., an additive process) X^ K In the 
time-independent case we have that (j)n{C} — 4'n{t, ^) is the characteristic exponent of the Levy process X'^"-^ 



Example 5 (Taylor series expansion). iPagliarani. Pascucci. and Rigal (|2013[) approximate the drift and 
diffusion coefficients of ^l as a power series about an arbitrary point a; € M. In our more general setting, 
this corresponds to setting Bn{x) ~ {x — x)" and expanding the diffusion and killing coefficients a{t, •) and 
7(t, •), as well as the Levy measure i'{t, ■, dz) as follows: 



a(i, x) = ^ an{t, x)Bn{x), 

oo 

7(^,2:) = ^7„(i,x)B„(x), 



n=0 

00 



an{t,x) = — 9"a(t, x), 



ln{t,x) = —d^-f{t,x), 



(18) 



i'{t,x,dz)=} i^nit,x,dz)Bn(x), Vn{t,x,dz) = —;d"i'(t,x,dz). 
^ — ' n\ 

In this case, (|T2|) and ([16]) become (respectively) 

CO CO 

(/)(i,a:,C)=^(x-S)>„(i,0, yi = ^(a;-.T)>„(i,2)), 

n=0 n=0 

where, for all n > 0, the symbol 0„(i,^) is given by P^ with coefficients given by ^TE\\ . The choice of 
X is somewhat arbitrary. However, a convenient choice that seems to work well in most applications is to 
choose X near Xt, the current level of X. Hereafter, to simplify notation, when discussing implementation 
of the Taylor-series expansion, we suppress the x-dependence: a„(t, x) — ?> a„(t), 7„(i,a;) — > 7„(t) and 

Vn{t,X,dz) — > Vn{t,dz). 

Example 6 (Two-point Taylor series expansion). Suppose / is an analytic function with domain M and 
a;i,a;2 G K. Then the two-point Taylor series of / is given by 



n=0 



where 



Co{xi,X2) 



f{x) = ^ {Cn{xi,X2){x- Xi) +Cn{x2,Xi){x - X2)) {x- Xi)"'{x- X2) 



(19) 



X2 — X\ 



,. _. ^ (fc + n - 1)! {-lYkd--'f{x,) + {-ir+^nd-;'f{x2) ,_ 
Cn(xi,X2) = 2^7:077 —, 7^ m+TTT^ • (20) 



fe=0 



k\n\{n — k)\ 



(Xl -X2)''+"+l 



For the derivation of this result we refer the reader to 



Estes and Lancasteij (|l966r ): lLopez and Temmd (|2002i ) 



Note truncating the two-point Taylor series expansion (llOp at n = m results in an expansion which of / 
which is of order ©(x^'"^^). 

The advantage of using a two-point Taylor series is that, by considering the first n derivatives of a 
function / at two points xi and X2 , one can achieve a more accurate approximation of / over a wider range 
of values than if one were to approximate / using 2n derivatives at a single point (i.e., the usual Taylor series 
approximation) . 

If we associate expansion (|19p with an expansion of the form /(x) — '^n=o fn^nix) then foBo{x) = 
c„(xi,X2)(x — xi) -|- c„(x2,xi)(x — X2), which is ajjine in x. Thus, the terms in the two-point Taylor series 
expansion would not be a suitable basis in (|12p since Bq{x) 7^ 1. However, one can always introduce a 
constant M and define a function 



^(x) := /(x) - M, 



so that 



!{x) = M + F{x). 



(21) 



Then, one can express / as 

oo 

fix) =M + J2 {CnM^i,X2){x - xi) + C„_i(a;2, Si)(x - x^)) {x - x^^x - a;2)"-\ (22) 



where the C„ are as given in I^U^ with f ^f F. If we associate expansion ([2^ with an expansion of the form 
f{x) — X]n=o fnBn{x), then we see that foBo{x) = M and one can choose ^0(2;) = 1- Thus, as written in 
the terms of the two-point Taylor series can be used as a suitable basis in (|12p . 
Consider the following case: suppose a{t,x), "/(t,x) and v{t,x,dz) are of the form 



iy{t,x,dz) — f{x)J^{t,dz), 



a{t,x)^f{x)A{t), j{t,x)^f{x)T{t), 

so that (t){t,x,C) = f{x)^{t,C) with 

N(t, dz)(e^ - 1 - z)iS^ + f J^it, dz){e'^^ - 1 - i^z) 



(23) 



It is certainly plausible that the symbol of A would have such a form since, from a modeling perspective, it 
makes sense that default intensit y, volatility and jump-inte n sity would be proportion al. Indeed, the Jump-to- 
default CEV model (JDCEV) of ICarr and Linetskvl (|2006r ): ICarr and MadanI (J2010l ) has a similar restriction 
on the form of the drift, volatility and killing coefficients. 

Now, under the dynamics of ([^5]) . observe that (j){t,x,£,) and A can be written as in p^ and ^ 
respectively with Bq — 1 and 



n > 1. 



(24) 



Bn{x) = {Cn-l{xi,X2){x- Xi) + C„-i{x2,Xi){x- X2))ix -Xi)"^ (x - ^2)" , 

As above C„ (capital "C") are given by ((20| with f ^f F := f — M and 

Mt, - A/$(t, 0, Mt, = Ht, 0, n>l. 

As in example [5l the choice of xi, X2 and M is somewhat arbitrary. But, a choice that seems to work well 
is to set xi — Xt — A and X2 ~ Xt + A where A > is a constant and M — f{Xt). It is also a good 
idea to check that, for a given choice of xi and X2, the two-point Taylor series expansion provides a good 
approximation of / in the region of interest. 

Note we assumed the form ((221) only for sake of simplicity. Indeed, the general case can be accommodated 
by suitably extending expansion (|12p to the more general form 

oo 3 

(l>{t,X,£,) = ^^Bi^n{x)(t)i,n{t,£.)i 
n=0 i=l 

where 4'i^n for z = 1,2,3 are related to the diffusion, jump and default symbols respectively. For brevity, 
however, we omit the details of the general case. 

Example 7 (Non-local approximation in wheighted L^-spaces). If {Bn}n>o is a fixed orthonormal basis in 
some (possibly weighted) L'^{W,m{x)dx) space, then in the expansion P^ we have 

^„(i,O = (0(t,-,O,^n(-))m, n>0. 
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A typical example would be to choose Hermite polynomials iJ„ centered at x as basis functions, which (as 
normalized below) are orthonormal under a Gaussian weighting 

B^ix)^H^ix-x), g„(x) :^ ^^^ ^""^P^7'\ n>0. 

In this case, we have 

Mt,0 ^ {(j)it,-,^),Bn)m ■= / (l3{t,x,^)Bnix)m{x)dx, m(x) := cxp ( - (x - x)^) , 

Once again, the choice of x is arbitrary. But, it is logical to choose x near Xt, the present level of the 
underlying X. In Figure [T] we compare the 2nd order two-point Taylor series, the 4th order (usual) Taylor 
series and the 4th order Hermite polynomial expansion for an exponential function. 

Now, returning to Cauchy problem ©, we suppose that v = v{t, x) can be written as follows 

oo 



oo 
V = 

n=0 



Following IPagliarani et al.l (|2013f ), we insert expansions ([T6| and p5|) into Cauchy problem ([6]) and find 

{dt + Ao)vo = 0, vo{T, x) = h{x), (26) 

n 
{dt+Ao)Vn^-^Bk{x)AkVn-k, Vn{T,x)^0. (27) 

fc=l 
We are now in a position to find the explicit expression for w„, the Fourier transform of v„ in (|26p - (|27p . 

Theorem 8. Suppose h G L^(R,dx) and let h denote its Fourier transform. Suppose further than u„ and 
its Fourier transform «„ exist. Then Vn{t,^) is given by 

voit,^) = exp ( / 0o(s,C)ds 1 m, (28) 



Vn{t,^)=^ CXp I / (j)o{u,l 



Odu]Bkiid^)(j)kis,Ovn~kis,Ods, n > 1. (29) 

Note that the operator Bk{id^) acts on everything to the right of it. 

Proof. See Appendix 1X1 D 

Remark 9. Assuming w„ G L^{M.,dx), one recovers u„ using 

vjt,x)^ [ d^-^e'i-v,,{t,(,). (30) 

JR V 27r 

As previously mentioned, to obtain the FK transition densities p{t,x;T,y) one simply sets h{x) = Sy{x). In 
this case, /i(^) becomes V'y(— 0- 

When the coefficients (a, 7, v) are time-homogeneous, then the results of Theorem [8] simplify considerably, 
as we show in the following corollary. 
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Corollary 10 (Time-homogeneous case). Suppose that X has time-homogeneous dynam,ics with the local 
variance, default intensity and Levy measure given by a{x), ^{x) and ^{x^dz) respectively. Then the symbol 
4>n{t,0 — 4'n{S.) is independent oft. Define 

T{t) —T -t. 

Then, for n <Q we have 

Vn{t,x) = Un{T(t),x) 

where 

So(r,0 = e-^"'^'^(0, (31) 

n .r 

Proof. The proof is an algebraic computation. For brevity, we omit the details. D 

Example 11. Consider the Taylor density expansion of Example El That is, Bn{x) = {x — x)". Then, 
explicitly, ui{t,^) and U2{t,^) are given by 

ui(t,e) = e*^"(«' (thiOxMO + ^tMOh'iO + l^t'HOMOMO + mOMo) , (32) 

- 2th'{0<f>2iO - t'mc^oiO<f>2{0 ~ \f^i{if^h!\i) ~ tMOh"{0 - It'hiO'PiiOKiO 

Higher order terms are quite long. However, they can be computed quickly and explicitly using the Mathe- 
matica code provided in Appendix IbI 

Remark 12. Many common payoff functions (e.g. calls and puts) are not integrable: h ^ L^(R,dx). Such 
payoffs may sometimes be accommodated using generalized Fourier transforms. Assume 

h{£_) := / dx e^^^^h{x) < oo, for some S, = £_r + i£.t with ^r,£.i G K. 

Jm v27r 

Assume also that (j){t,x,£^r + iS,i) is analytic as a function of ^r- Then the formulas appearing in Theorem 
[8] and Corollary [TOl are valid and integration in (l30t is with respect to S,r (i-c, dS, -^ rf^r)- For example, the 
payoff of a European call option with payoff function h{x) — (e^ — 6*"')+ has a generalized Fourier transform 

M0= /rf^-^e^'«-(e--e'=)+ = -=^^-^, ^ = Cr + ^C., CreR, C,e{-^,-l). 
Jm ^^"^ V27r(i^ + ^^) 
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In any practical scenario, one can only compute a finite number of terms in (|25|) . Thus, we define v^^', 
the A'^th order approximation of v by 

N N 

The function u^^' {t, x) (which we use for time-homogeneous cases) and the approximate FK transition 
density p^^'{t, x; T, y) are defined in an analogous fashion. 



Remark 13. As in 



Pagliarani et al. 



()2013l) . when considering models with Gaussian-type jumps, i.e., models 
with a state-dependent Levy measure v(t, x, dz) of the form (1341) below, all terms in the expansion for 
the transition density become explicit. Likewise, for models with Gaussian-type jumps, all terms in the 
expansion for the price of an option are also explicit, assuming the payoff is integrable against Gaussian 
functions. Indeed, as we shall see in Appendix [Cl the main term p*-"' {t, x; T, y) is given by a series of Gaussian 
densities; further approximations p^^' {t, x;T,y) can be represented as a suitable differential operator applied 
to p^^\t,x]T,y) acting on the variable x (see for instance Proposition l311 for the first order). 

5 Pointwise error bounds for Gaussian models 

In this section we prove some pointwise error estimates for p^^'{t, x; T, y), the A^th order approximation of 
the FK density of {dt +A) with A as in (|S]). Throughout this Section, we assume Gaussian-type jumps with 
(t, a;)-dependent mean, variance and jump intensities. Furthermore, we work specifically with the Taylor 
series expansion of Example [S] That is, we use basis functions Bn{x) — {x — x)" . 

Theorem 14. Assume that 

m<a{t,x)<M, 0<-f{t,x)<M, t e [0,T], a; G M, 

for some positive constants m and M , and that 

\(t X) (»-Ai(t,x))2 

i^{t,x,dz) = A(t,a;):N-^(t,,),52(t,,)((iz) := ' ' e =^"('.-) dz, (34) 

\'2'Kd[t,x) 



with 

m<S'\t,x)<M, < A(i,x),|/i(i,x)| < M, t e [0,T], x eM.. 

Moreover assume that a,^, X,6, fi and their x-derivatives are bounded and Lipschitz continuous in x, and 
uniformly hounded with respect to t ^ [Oi^]- Let x ~ y in psp . Then, for N > 1, we hav^ 

Pit, x; T, y) - p^''\t, x; T,y)\< g^iT - t) (f (t, x; T, y) + ||9,i/||^ f{t, x; T, y)) , (35) 

for any x, y e M and < t < T, where 

gN{s) = (s) , as s -^ 0+. 

^Here ||9a:2^||oo = '^^'^i 1 1 ^^a;-^ I loo ' ll'^^'^llcxD ' ll'^^A'llcx)}' where || ■ ||oo denotes the sup- norm on [0,T] X R. Note that ll^a;!^!!;,^ = 
if A, 5, fi are constants. 
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Here, the Junction T is the fundamental solution of the constant coefficients jump-diffusion operator 

M - f 

dtu{t,x) + —dxx + M / {u{t,x + z) - u{t,x))>iM^M{dz), 



where M is a suitably large constant, and T is defined as 

fit, X- T,y) = f2 M!^!(^_^e'=+^f (i, x; T, y), 

' ^ A/ft*' 

k=Q ^ft- 

and where C being the convolution operator acting as 

e/(x)= f f{x + z)J^a.Midz). 
Jr 

Proof. See Appendix [Cl D 

Remark 15. In the proof of Theorem [T4l we see that functions C'^F take the following form 

e^f (^ x-Ty)- e-^^(-*) f ^ (^^(^-^))" exp ( J- - V + ^^n + k))' ] ^ ^ ^ 

^' ' '^^" t',nW2nM{T-t + n + k) ""y 2M{T - t + n + k) ) ' " '' 

and therefore F can be explicitly written as 



r{t,x;T,y) = e-^^('-'' > , ,. , ' '_ " . ^ exp 



(M(r-i))"^^ / (a;-y + M(n + fc + l))' 



.A:=0 



!V^V27rA/(r - t + n + A: + 1) V 2Af (T - t + n + A: + 1) 



By Remark [T5l it follows that, for any x ^ y fixed, we have T{t, x; T, y) = 0{T — t) as (T — t) tends to 
0, whereas G'^T{t,x;T,y), and thus also r{t,x;T,y), tends to a positive constant as (T — t) goes to 0. It is 
then clear by ([35]) that, with x ^ y fixed, the asymptotic behavior of the error, when t tends to T, changes 
from (T — t) to [T — t)'^ depending on whether the Levy measure is lo cally-dependent or not. 



Theorem [T4l extends the previous results in iPagliarani et alJ (|2013[ ) where only the purely diffusive case 



(i.e A = 0) is considered. In that case an estimate analogous to ([55| holds with 

gAr(s) = (s^5~ j , ass^O+. 

Theorem [T3] shows that for jump processes, increasing the order of the expansion for N greater than one, 
theoretically does not give any gain in the rate of convergence of the asymptotic expansion as i — >■ T~ ; this 
is due to the fact that the expansion is based on a local (Tayl or) approximation while the FIDE contains a 



non-local part. This estimate is in accord with the results in iBenhamou et al.l (|2009l ) where only the case 
of constant Levy measure is considered. Thus Theorem [T3] extends the latter results to state dependent 
Gaussian jumps using a completely different technique. Extensive numerical tests showed that the first 
order approximation gives extremely accurate results and the precision seems to be further improved by 
considering higher order approximations. For example, in Figure [5] we plot the approximate transition 
density p^^\t, x; T, y) for different values of N for the Levy-type model considered in Section [7. 1.11 We note 
that, for T — t < 5.0 years, the transition densities p^^'^ (t, x; T, y) and p'^'^^ {t, x; T, y) are nearly identical. This 
is typical in our numerical experiments. 

Remark 16. The proof of Theorem [14] is also interesting for theoretical purposes. Indeed, it actually 
represents a procedure to construct p{t,x;T,y). Note that with p^'^'{t,x;T,y) being known explicitly. 



Equation ([35]) provides pointwise estimations for the exact FK density as well. 
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6 Implied volatility for local volatility models 

For European calls and puts, it is the implied volatility induced by an option price - rather than the option 
price itself- that is the quantity of primary concern. In this section we derive an implied volatility expansion 
for time- homogeneous local volatility models (i.e., for models with no jumps and no possibility of default). 

Assumption 17. In this section only, we assume 



St = e^\ dXt = -a{Xt)dt + y^2a{Xt)dWt, 

so that 00 (C) = ao(-C^ - ^0 in dUl), with oq > 0. 

Remark 18. Note that in the Taylor expansion of Example [SJ we have oq = a{x). In the two-point 
Taylor expansion of Example |6l we have ao — MA. And in the Hermite expansion of Example [7l we have 

oo = /r a(2;)e^'-^^^-' dx. 

Throughout this section, we fix an initial value of the underlying Xq — x, a time to maturity t and 
a call option payoff h{Xt) — (e^* — e'^)+. The price of this option u{t,x) — u{t,x;k) can be computed 
(approximately) using Corollary 1101 To simplify notation, throughout this section we will suppress all 
dependence on (i, x, k). However, the reader should keep in mind that option prices and their corresponding 
implied volatilities do depend on these variables. We begin our analysis with a definition of the Black-Scholes 
price. 

Definition 19. The Black-Scholes Price u^^ : R+ -^ K+, defined as a function of volatility cr, is given by 

Note that 4>^ {-'jCr) is simply the Levy exponent of a Brownian motion with volatility a and drift ~^cr^. 
Thus, this is simply the Fourier representation of the usual Black-Scholes price. 

Remark 20. Note that, by equations (pO)) and (|3T| . under Assumption [iTl we have uq = u^^W'^o.o)- 

Definition 21. The implied volatility induced by option price u is the unique positive solution a of the 
equation 

u^^(cr) = u. (36) 

We are now in position to derive an asymptotic expansion for implied volatility. We begin by writing the 
option price and implied volatility as follows 



oo 
U — 

n=0 



X]e"«n, e = l, (37) 



n=l 
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We introduce the constant e — 1 purely for the purposes of accounting; e plays no role in our final result. 
Next, we expand u^^ as a power series about the point ctq. We have 



,BS 



oo oo / oo \ 



00 ^ 

^^"(-0) + E ^ 

n— 1 
00 

.^^(ao)+E=' 



fe=i 



E E n-.. 

k—n \ji-\ hjn— fc i— 1 



E;^ E n^.h^ 






,BS 



M + E=' 



fe=l 



n. , 

n— 1 \jiH {-jiz—ki—l 

"=2 \jl+---+i„=fci=l 



zi"^(ao). 



(38) 



Now, we insert expansions (|57| and ([55)1 into (pS)) and collect terms of like order in e 
0(1): uo^u^'^iao), 



0(£'=): ^,^a,5.«5^(ao) + E;^ E n'^.0^>''^( 

n=2 ■ yji+...+j„=fc4=l 

Solving the above equations for {crfe}^o ^*^ ^'^'^ 
0(1) : ao = V2^, 

k 

0(e*= 



Co) 



1 



CTfc 



9.ii^^(co) 



"-eM e n-.>>"'M 



fc> 1. 

(39) 
fc > 1, (40) 



n=2 \ji+...+j„=A;i=l 

where we have used uq — u^^{y/2ao) to deduce that ctq = V^ao (see Remark [^ . Observe, the right hand 
side of POJ) involves only cr^- for j < k — 1. Thus, the {(JkJkLi can be found recursively. Explicitly, up to 
O(e^) we have 

Ml 



0(e): C7i = 

daUQ 



1 2 n? 
■"2 - 2TO'iC'^^«0 



0(e3): as 
O(e^) : a4 



9crU0 

d„Uo 

U4 - (asgi^g + icrfgg + ^g2CTi9g + ^gfa;^)Mo 

dcrUo 



(41) 



where 



Theoretically, one can use recursion relation (j40p ad infinitum to find an for arbitrarily large n. To sum up, 
our nth order asymptotic approximation for implied volatility is as follows: 

n 
fc=0 
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where the {cfc} are given by ([M|) and (|40l) . 

As written, equations (j4ip are not very convenient. Indeed, for each u„ and for all terms of the form 
d^uo, one must be compute Fourier integrals. However, with a little more work, we can compute the {cr,} 
explicitly with no integrals and no special functions. The first step is to express u„ as 

Un^Luo, £ = ^bk{t,x)d^{dl-dx), 

fe>0 

where the {bk} are functions of {t,x). 

Assumption 22. For the remainder of Section [51 we assume Bn{x) and 0„(^) = (t>n{tX) £^s in Example [S] 
That is, Bn{x) = [x — x)" and cf)n{C} = an(^C^ ^ *C) with a„ = ^c)"a(i). This choice is not necessary, but 
simplifies the computations that follow. Similar implied volatility results can be derived for other choices of 
Bn{x) and 0„(O. 

As shown in Appendix ID| under Assumption 1221 for ui we have 

ui{t,x) = {t{x-x) + \fao{2d^-l)Yi{dl~d.,)uQ{t,x). (42) 

and for U2 we have 

U2{t,x) = it{x - xf +f{x- x)aQ{2d^ - 1) + lt^al{2d^ - 1)^ + i^ao jazia^ - d^)uQ{t, x) 

+ (\t\x - xfa^{dl - d,) + ii^ (^ _ x) a^idl - d,)ao{2d, - 1) 



+ ^t^aiidl - d,)ali2d, - If + le {x - x) a^{2d., - 1) 

+ \t\o{2d^ - l)ai(29, - 1) + \t\^{dl - d^)a^ai{dl ~ d^)uo{t,x). (43) 

Now, using equations P^ - P5)) . as well as 

9™(a2 - d^)u^'^{a) __ ^"exp (x - %) 



d^u^'^ia) = ta{dl - d^)uB%a) 



(9^ - a.)«^^(a) exp(x-f) 



where (i_|_ = — ^ (a; — fc + 2icr^), a straightforward (but tedious) algebraic computation shows that 



(72 == - 



daUQ 
02 



f fc^ + /sx + x^ + too - 3(fc + a;)x + 35^] 

2 

°^^ (too(6 + too) + 6 (fc2 + x^ - 2(A: + x)x + 2x^) \ . 



3\/2ao 

48ao\/2ao 
If we set a; — >■ a; then we have the following second order approximation for implied volatility 



a'-'^' = (To + (Tl + 0-2 

I o-i 

= V 2an 



2 

-.{k-x) + -^^({k-xf+tan) ~ °^ U{k - xf + too (6 + too) ) , (44) 



2\J2aQ 3v2ao V / 48ao\/2ao 

which is quadratic in (A; — x) and in t. 
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7 Examples 

In this section, in order to illustrate the versatility of our asymptotic expansion, we apply our approximation 
technique to a variety of different Levy-type models. We consider both finite and infinite activity Levy- 
type measures, models with and without a diffusion component, and models with and without jumps. We 
study not only option prices, but also implied volatilities and credit spreads. In each setting, if the exact or 
approximate option price/implied volatility/credit spread has been computed by a method other than our 
own, we compare this to the option price/implied volatility/credit spread obtained by our approximation. 
For cases where the exact or approximate density/option price is not analytically available, we use Monte 
Carlo methods to verify the accuracy of our method. 

Note that, some of the examples considered below do not satisfy the conditions listed in Section [2l In 
particular, we will consider coefficients (a, 7, v) that are not bounded. Nevertheless, the formal results of 
Section m work well in the examples considered. 

7.1 CEV-like Levy- type processes 

We consider a Levy- type process of the form ([l} with CEV-like volatility and jump-intensity. Specifically, 
the log-price dynamics are given by 

a(x) = i(52g2(^-i).^ z.(a;,dz)=e'('3-i)"?^(dz), lix) = 0, S > 0, /3e[0,l], (45) 



where J^^dx) is a Levy measure. When 3\f = 0, this model reduces to the CEV model of lCoxl (|1975l ). Note 
that, with j3 g [0,1), the volatility and jump-intensity increase as x -^ — cxi, consistent with the leverage 
effect (i.e., a decrease in the value of the underlying is often accompanied by an increase in volatility /jump 
intensity). This characterization will yield a negative skew in the induced implied volatility surface. As the 
class of models described by (|45|) is of the form (|23p with f{x) — e^^^~^^^, this class naturally lends itself 
to the two-point Taylor series approximation of Example [S] Thus, for certain numerical examples in this 
Section, we use basis functions _B„ given by (j24p . In this case we choose expansion points xi and X2 in a 
symmetric interval around Xq and in (I^Tt we choose M = ,f{Xa) = 3"^^^"^^-^°. For other numerical examples, 
we use the (usual) one-point Taylor series expansion B„{x) — {x — 5)". In this cases, we choose x = Xq. 
We will consider two different characterizations of ?sf(dz): 



Gaussian: J'i{dz) = A^== exp ( ^ ^ — — ) dz, (46) 



-(z — m~*^ 
V27rs2 "" V 2s- 

g-A_|2| g-A 

Variance- Gamma: l^{dz) = ( - — — j — l{z<o} ^ I{z>o} ) d^^ (47) 




A.= h/^ + — ±-j 

Note that the Gaussian measure is an example o f a fini te-activity Levy measure (i.e., 3\f(R) < 00), whereas the 
Variance-Gamma measure, due to lMadan et al.l (|l998l ). is an infinite-activity Levy measure (i.e., N(R) = 00). 



As far as the authors of this paper are aware, there is no closed-form expression for option prices (or the 
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transition density) in the setting of (|^5|) . regardless of the choice of 'N{dz). As such, we wih compare our 
pricing approximation to prices of options computed via standard Monte Carlo methods. 

7.1.1 Gaussian Levy Measure 

In our first numerical experiment, we consider the case of Gaussian jumps. That is, J^{dz) is given by (j46p . 
We fix the following parameters 

,5 = 0.20, /3 = 0.25, A = 0.3, to = -0.1, s = 0.4, 5o = e"^ = 1. 

Using Corollary[TUl we compute the approximate prices m'^'^' {t, x; K) and w'-'^-' (t, x; K) of a series of European 
puts over a range of strikes K and with times to maturity t = {0.25, 1.00, 3.00, 5.00} (we add the parameter 
K to the arguments of u'"-* to emphasize the dependence of u*^"-' on the strike price K). To compute 
u^*-'(i,x; iC), i — {0,3} we use the we the usual one-point Taylor series expansion (Example [5]). We also 
compute the price u{t, x; K) of each put by Monte Carlo simulation. For the Monte Carlo simulation, we use 
a standard Euler scheme with a time-step of 10~^ years, and we simulate 10^ sample paths. We denote by 
y(A^c) ^^^ ^. j^^ ^j^g price of a put obtained by Monte Carlo simulation. As prices are often quoted in implied 
volatilities, we convert prices to implied volatilities by inverting the Black-Scholes formula numerically. That 
is, for a given put price u{t, x; K), we find a{t, K) such that 

u{t,x]K) ^u^'^{t,x]K,(T{t,K)), (48) 

where u^^{t, x; K, a) is the Black-Scholes price of the put as computed assuming a Black-Scholes volatility 
of (7. For convenience, we introduce the notation 

IY[u(t,x;K)] :^a{t,K) 

to indicate the implied volatility induced by option price u{t, x; K). The results of our numerical experiments 
are plotted in Figure [31 We observe that l\[u^^\t, x; K)] agrees almost exactly with TV[u''^'-'\t, x; K)]. The 
computed prices u^-^^t, x; K) and their induced implied volatilities lV[u^^^t, x; K)], as well as 95% confidence 
intervals resulting from the Monte Carlo simulations can be found in Table [TJ 

7.1.2 Variance Gamma Levy Measure 

In our second numerical experiment, we consider the case of Variance Gamma jumps. That is, 1^{dz) given 
by (|47p. We fix the following parameters: 

5 = 0.0, /3 = 0.25, 6' = -0.3, p = 0.3, k = 0.15, 5o = e^ = 1. 

Note that, by letting 5 — Q, we have set the diffusion component of X to zero: a{x) ~ 0. Thus, X is 
a pure-jump Levy-type process. Using CoroUarv llOi we compute the approximate prices u^ \t,x;K) and 
u'^'^\t^x]K) of a series of European puts over a range of strikes and with maturities t G {0.5,1.0}. To 
compute u^^\ i e {0,2}, we use the two-point Taylor series expansion (Example IH]) • We also compute the 
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put prices by Monte Carlo simulation. For the Monte Carlo simulation, we use a time-step of 10 "^ years 
and we simulate 10^ sample paths. At each time-step, we update X using the following algorithm 



Xt+At ^Xt + b{Xt)At + j+{Xt) - -/-(Xt), 



b{x) 



I{x) 



loa 



A_ 



1 + A_ 



log 



A, 



A+-1 



7±(a;)^r(/(a;).At/K,l/A±), 



where T{a, b) is a Gamma-distributed random variable with shape parameter a and scale parameter b. Note 
that this is equivalent to considering a VG-type process with state-dependent parameters 



k'(x) :— k/I{x), 



9'{x) :=6'/(x), 



p'ix):^py/l{^. 



These state-dependent parameters result in staie- independent X± (i.e., A± remain constant). Once again, 
since implied volatilities rather than prices are the quantity of primary interest, we convert prices to implied 
volatilities by inverting the Black-Scholes formula numerically. The results are plotted in Figure ID We 
observe that TV[u^^^t,x;K)] agrees almost exactly with l'V[u'^^'^^\t,x;K)]. Values for u'^'^\t,x]K), the 
associated implied volatilities lY[u^'^\t,x]K)\ and the 95% confidence intervals resulting from the Monte 
Carlo simulation can be found in table [2l 



7.2 Comparison with iLorig ( l2012l ) 

In 



Lorisj (|2012l ). the author considers a class of time- homogeneous Levy- type processes of the form: 

a{x)^-{bl + eblr^{x)), 1 



7(x) = Co + eci?7(a;), 
v{x, dz) — iyo{dz) + eri{x)h'i{dz). ^ 



(49) 



Here, {bo,bi,co,ci,e) are non-negative constants, the function ry > is smooth and vq and vi are Levy 
measures. When ri{x) = eplx) := e^^, the author obtains the following expression for European-style 
options written on X 



u{t,x) = y^£"w„(^,x), 

n{t,x) = enfi{x) I C?^ X! fF 



(50) 



gt'Ji'e-ifc^ 



fn-l 



where t is the time to maturity, x is the present value of X and 



n^«-»'=/3 ]MoMx)- 



\k=Q 



T^i - ^^0 {-(" - *C) + cM -I)- I Mdz)(e'' - 1 - z)ze + / Mdz)(e'^' - 1 - i^z) , 

1 {dz) fe^ - 1 - zW + I vx {dz) ('e*^^ - 1 - i^z) ■ 



As in Corollarv [TUl h{^) is the (possibly generalized) Fourier transform of the option payoff h{x), and -0^ is 
as given in ([TI 
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Now consider the following model: 



a{x)^Af{x), 



lix) = r/(a;), 



v{x,dz) = f{x)J<idz) 



The models described by P^ and ([51]) coincide if we choose 
1, 



-bf = a,A, 



Ci — GiT, 



Vi{dz) ~ Qi'Nldz), 



f{x) ^ ao + eair]{x), (51) 



z = {0,l}. 



Furthermore, comparing equations (j23p with (|5ip . we see that (j5ip is precisely the form considered in 
Example m Thus, in this Section we use the two-point Taylor series approximation of Example |6] with basis 
functions i?„ given by (1241) . We choose expansion points xi and X2 in a symmetric interval around Xq and 
in dUl) we choose M = /(Xq) = e/^-^o. 

In our numerical experiment, we consider Gaussian jumps (i.e., J^{dz) given by (j46])) and we fix the 
following parameters: 



s — ao — ai ^ 1, 
m = -0.2, 



/3 = -2.0, 
t = 0.5, 



A^-0.15^ 
2 

Xo = 0.0, 



r = o.o, 

xi = —0.3, 



A = 5 = 0.2, 
X2 = 0.3. 



Using Corollary[TUl we compute the approximate prices m'^'^^ (i, x; ii') and u^^' (i, x; ii') of a series of European 
puts with strike prices K S [0.5, 1.5] (once again, we add the parameter K to the arguments of u^") to 
emphasize the dependence of u'"' on the strike price K). We also compute the price u(i, x; K) using ([50l) . In 
()50p . we truncate the infinite sum at n = 8. As in Section F7.H we convert option prices to implied volatilities. 
The results are plotted in Figure [51 We observe a nearly exact match between the induced implied volatilities 
IV[u^^^(i, X] K)] and IV[u(i, x; K)], where u(i, cc; K) (with no superscript) is computed by truncating ([50)) at 
n = 8. 



7.3 Comparison to NIG-type processes 



A Normal Inverse Gaussian (NIG) (see iBarndorfii-Nielsenl ( 19981) ) is a Levy process with characteristic Levy 
exponent 4'{^) given by 



4>{0=^^li-5 ^o? - (/3 + liY - ^M^^ 



In Chapter 14, equation 14.1 of iBovarchenko and Levendorskiil (|2000[ ). that authors consider NIG-like Feller 
processes with symbol 



<i){x,i)=iii{x)i-8{x) ^Ja^(x) - (/3(x) + i^Y - yJo?(x) - /32(a;) 
where /i,(5, a, /3 e C/J"(M), (5, a > 0, /^, /3 e M, and where there exist constants c and C such that 5(x) > c. 



a(a:) — |/3(x)| > c and |/i(2;)| < C Note that if X is a NIG-type process with symbol ^{x,C)i then S 



^,x 



is a martingale if and only if (j){x. 



0. Thus, the triple {a, (3,5) fixes ^. 



Bovarchenko and Levendorskiil ([20001 ) deduce the following asymptotic expansion for u{t, x) (see the 

m. (52) 



equations following (14.27) and equation (16.40)). 

u{t,x) ■.^¥.MXt) = / dc4=e^«"e**(-^«) (l + ]-e[zd^^{x,Om<l^{^M + 
JR v27r V ^ 
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We note that, if one chooses basis functions Bn{x) — (x— x)" as in Example[S]with x — x, then (/)(x, ^) = (/)o(C) 
and dx(f>{x,^)d^<j){^TO — 0i(C)'/'o(C)- Thus, from Corollary [TOl and equations (|3T|) and ([32|) . it is easy to see 
that expansion (j52p is contained within the first two terms of the (one-point) Taylor expansion obtained in 
this paper. 



7.4 Yields and credit spreads in the JDCEV setting 

Consider a defautable bond, written on S, that pays one dollar at time T > t if no default occurs prior to 
maturity (i.e., St > 0, C > T) and pays zero dollars otherwise. Then the time t value of the bond is given 

by 



Vt - E[I{^>T}|yt] = Iic>t}vit,Xt;T), 



v{t,Xt;T) ^E[e- 1^^ ^'^'■^^''>'^'\Xt 



We add the parameter T to the arguments of v to indicate dependence of v on the maturity date T. Note 
that v(t,x;T) is both the price of a bond and the conditional survival probability: Q(C > T\Xt — x,(^ > t). 
The yield Y{t, x; T) of such a bond, on the set {C, > t}, is defined as 

-logvit,x;T) 



Y{t,x;T):= 



T-t 



(53) 



The credit spread is defined as the yield minus the risk- free rate of interest. Obviously, in the case of zero 
intere st rates, we have: y i eld = credit spread. 



In 



Carr and Linetskvl (|2006r ). the authors introduce a class of unified credit-equity models known as 



Jump to Default Constant Elasticity of Variance or JDCEV. Specifically, in the time-homogeneous case, the 
underlying S is described by ([1]) with 



a{x) = ^S^e'f^^, 



^{x) — b + cS e 



2 2Px 



v{x, dz) = 0, 



where 5 > 0, 6 > 0, c > 0. We will restrict our attention to cases in which /3 < 0. From a financial perspective, 
this restriction makes sense, as it results in volatility and default intensity increasing as S* — > 0+, which is 
consistent with the leverage effect. Note that when c > 0, the asset S may only go to zero via a jump from a 
stric tly positive value. That is, ac cording to the Feller boundary classification for one-dimensional diffusions 
fsee lBorodin and SalminenI (|2002| ). p. 14), the endpoint — oo is a natural boundary for the killed diffusion X 
(i.e., the probability t hat X reaches — oo in finite ti me is zero). The survival probability v{t,x\T) in this 

poiol ). equation (8.13). We have 



setting is computed in 



Mendoza-Arriaga et al. 



v(i, X] T) = uiT — i, x) 



E 



,-(ft+c.„)(T-t) r(l + c/|/3|)r(n + l/(2|/3|)) 



r(^ + l)r(l/(2|/3|))n! 
X Ai/(2|/5|)e^exp (-Ae-^/S-) iFi(l -n + cl\l3\;v + l-Ae-''^^)\ 

where iFi is the Kummer confluent hypergeometric function, r(a;) is a Gamma function and 
1 + 2c , b 



(54) 



A = 



6^ 



uj^2\B\b. 
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We compute u(T — t, x) using both equation (|54p (truncating the infinite series at n = 70) as well as using 
Corollary [TOl We use basis functions from the Taylor series expansion of Example [5] Bn{x) ^ (x — x)". 
After computing bond prices, we then calculate the corresponding credit spreads using (|53p . Approximate 
spreads are denoted 



y'-fex;r).^ -'°'f';---^' 



T-t 

The survival probabilities are and the corresponding yields are plotted in Figure [SI Values for the yields 
from Figure |6] can also be found in Table [3] 

Remark 23. To compute survival probabilities v{t, x; T), one assumes a payoff function h{x) = 1. Note that 
the Fourier transform of a constant is simply a Dirac delta function: h{£,) = S{£^). Thus, when computing 
survival probabilities and/or credit spreads, no numerical integration is required. Rather, one simply uses 
the following identity 



to compute inverse Fourier transforms. From the above identity and equations p2p - p3p one easily obtains 

«o(i, a;) = £-("+*' -'^')*, 






2 ^ '^8 ^2 ^2 ' 

7.5 Implied Volatility Expansion for CEV 



In this Section we apply the implied volatility expansion of Section [S] to the CEV model of lCoxl (jl975() . The 
log-price dynamics are given by 

a(a;) = ifj^e^t^-i)^, y{x,dz)^Q, l{x)=Q, 5 > 0, /3e[0,l], 

In this setting the exact price of a call option with strike K = e^ can be expressed as follows: 

u{t, X- K) = e-Q{K, 2 + 2^, 2x) - e^^ (l - Q(2x, 3!^, 2a.)) , (55) 

^^ (M/2)"e-^/^r(^/2 + n,W2) 
Q(u;,.,M) = \^ (^ -, r(./2 + n) 

2g(2-/3)2; 
X 



(52(2 -/3)2i' 

2e(2-/3)fc 
.52(2 -/3)2f' 
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where T(a) and T(a, b) are the complete and incomplete Gamma functions respectively. For a given call 
option, the (almost) exact implied volatility a corresponding to the CEV model can be obtained by solving 
(H5 | numerically, with ujt , x\K) given by (I55t . 



Hagan and WoodwardI (|l999l ) derive the following implied volatility expansion for the CEV model 



„„„^ . r (wK2+«/^^ ^iw)!^^^^.v ^1,^.^^.,^ ,,,, 



/1-/3 I 24 V / / 24 /2(i-0) i ' •' 2 

Our second order implied volatility expansion can be computed using equation (j44p . We compare the 2nd 
order implied volatility expansion to the Hagan- Woodward expansion (|56p in Figures [7] and |51 Over the 
strikes and maturities tested, our implied volatility expansion performs favorably. 

8 Conclusion 

In this paper, we consider an asset whose risk-neutral dynamics are described by an exponential Levy-type 
martingale subject to default. This class includes nearly all non-negative Markov processes. In this very 
general setting, we provide a family of approximations - one for each choice of the basis functions (i.e. Taylor, 
two-point Taylor, L^ basis, etc.) - for (i) the transition density of the underlying (ii) European-style option 
prices and their sensitivities and (iii) defaultable bond prices and their credit spreads. For the transition 
densities, and thus for option and bond prices as well, we establish the accuracy of our asymptotic expansion. 
We also derive, for local volatility models, an expansion for the implied volatility of European calls/puts. 
Finally, we provide extensive numerical examples illustrating both the versatility and effectiveness of our 
methods. 

A Proof of Theorem [8] 

Proof. The formal adjoint of An is defined as the operator A\-^ such that 



{f,-An9) = {Aif,g), {u,v) = / u{x)v{x)dx, u,v € 



The operator A\ can be obtained formally by integrating by parts, which results in making the change 
dx — >■ —dx- That is 

yi,\-0„(t,-B), 'S = ~idx, 

with (/)„ as in ([T2|) . Note also that by (TTS]) and p7)) we have 

yiJ,e-"« = 0„(t, Oe-"«, i?„(a;)e-'^« = Bn{id{)e-'^^. (57) 

In the particular case of the expansion in Example [5l we have explicitly 

■^i - ln{t){-dx - 1) + an{t){dl + dx)+ [ J^n{t,dz){e^' ^l^z)dx+ I Vn{t, dz)(e-^^^ - 1 + zdx). 
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We Fourier transform equatfon (P7|) . First, we focus on the left-hand side 

= a,«(i,e) + Mt.o I dx^e-''^^vn{t,x) (by m) 

^ JR 

Next, for the right-hand side of ([27| . by ([57]) we have 

- V" / dx ( ^-_ gfc(a;) ) ■/lfc«„_fc(t, a:) = -y^ / dx ^^(^9^)^== I yifeW„_fc(i, x) 
j,^;^7e V V Stt / fc=i"'R V V Stt / 

= -V'Sfc(ia^) / dxLA|,— == w„_fc(i, a;) 

^-J^Bkiid^) (Mt,0 / dx^«„_fc(i,x) (by m) 
fc^i V JR V2tt J 

n 



fc=l 



Thus, we have the foUowing ODEs (in t) for Vn{t,^) 

{dt + Mt,0)yo{t,C) = o, vo{T,o^m, (58) 

n 

(at-f-^o(i,0)^n(^,0 = -E^fc(*^«)(<^'«(*'^)^"-fe(*'^)) «„(T,O=0. (59) 

fc=i 

One can easily verify (e.g., by substitution) that the solutions of (|58|) and (l59t are given by ((28)) and (|29)) 

respectively. 

D 

B Mathematica Code 

The following Mathematica code can be used to generate the u„(t,^) automatically: 

11 

u[n_,t.,A_]:=Exp[t0o[A]]^ 



/ Exp[-s</)o[A]]( V 
•^0 V,„=c 



k=l 

k! 



-^ ni!(k — m)! 



(-x)''-'°i'°D [0k[A]u[n - k, s. A], {A, m}] ds; 



u[O,t.,A_]=Exp[t0o[A]]h[A]; 
The function u{t,n,X) is now computed explicitly by typing {{[n, t, A] and pressing Shift+Enter. 



C Proof of Theorem [14 



Proof. For sake of simplicity we only prove the assertion when the default intensity and mean jump size 
are zero 7 = /i = 0, when the jump intensity and diffusion component are time-independent a{t,x) = a{x), 
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A(i, a;) = A(x) and when the standard deviation of the jumps is constant 5{t, x) = 5. Thus we consider the 
integro-differential operator 

Lu{t,x) = dtu{t,x) -\ —{dxx - dx)u(t,x) - A(x) (e^ - ij dxu{t,x) 



+ A(x) / {u{t^x + z) — u{t,x)^vs2[dz), 



with 



vs2 (dz) 



2ttS 



'2s^dz. 



Our idea is to use our expansion as a par ametr i x. Th at is, our expansion wih as the starting point of 
the classical iterative me thod introduced by Levil ([1907|) to construct the fundamental solution p{t,x;T,y) 



Pagliarani et al 



(|2013f ). we take as a parametrix our A'^-th order approximation 



of L. Specifically, as in 

p^^-'(i,x; T, y) with basis functions Bn = {x — S)" and with x = y. W e fir st prove the case JV = 1. 



By an a logy with the c lassical approach (see, for inst ance. iFriedmanI (119641) and 



Di Francesco and Pascucci 



(60) 



(|2005[ ). lPascuccil (|201ll ) for the pure diff^usive case, or lGarroni and Menaldil (|l992f ) for the integro-differential 
case), we have 

p{t,x;T,y)^p<^'\t,x-,T,y)+ f f p<^^Ht,x;s,OHs,e,T,y)d^ds, 

Jt Jv. 

where $ is determined by imposing the condition 

{) = Lp{t,x;T,y)^Lp'^^\t,x-T,y)+ f f Lp<^°\t,x; s,^Ms,e,T,y)d^ds - ^t,x;T,y). 

Jt Jr 

Equivalently, we have 

^t,x;T,y) = Lp^'\t,x;T,y)+ f f Lp^°\t,x;s,OHs,e,T,y)d^ds, 

Jt Jr 

and therefore by iteration 

oo 

<^{t,x;T,y) = ^Zn{t,x;T,y), 

where 

Zo{t,x;T,y) := Lp^^\t,x;T,y), 

Zn+i{t,x;T,y) :^ / / Lp^°\t,x]s,S)Zn[s,tT,y)d^ds. 
Jt Jr 



(61) 

(62) 
(63) 



The proof of Theorem [14] is based on several technical lemmas which provide pointwise bounds of each term 
Zn in ([CTjl . These bounds combined with formula (IMl) give the estimate of \p{t,x]T^y) — p'^^^(t, a;;T, j/)!. 
For any a,9 > Q and ^ > 0, consider the integro-differential operators 

L°'''''^u{t, x) = dtu{t, x) + ^{dxx - dx)u{t, x)~l (ei - l) dxu{t, x) + £ j {u{t, x + z) - u{t, x))i^g{dz), 

L ' ' u{t,x) ~ dtu{t,x) + —dxxu(t,x) -\- i {u{t,x -\- z) — u{t,x))i'e{dz). 
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The function r"'^'^(i, x; T, y) := r"^^'^(T -t,x~y) where 



n! 

n=0 



^ — ^ n' 



r^'«-^(i,x):=-==L===cxp 



is the fundamental solution of L"'^'^. Analogously, the function f "'^^^(i, x; T, y) := r^'^'^lT — t,x — y) where 

n=0 

f°-^(t,x):^-= /, ^ exp 



^27r(at + n6') V 2{at + nO) J ' 

is the fundamental solution of L ' ' . Note that under our assumptions, at order zero we have 

p^''\t,x;T,y)^T^^y^''''^<^y\t,x;T,y). (64) 

We also introduce the convolution operator Cg defined as 

Gefix) = f fix + z)^=e-idz. (65) 

Note that, for any 9 > 0, we have 



e,r-^'^(V)(x) = e-^*$:i^r:;Y(i 



X 
n=0 



e,f-«-^(t,.)(x) = e-"f;M:f:;^(i,x) 



n< 

In the rest of the section, we will always assume that 

m<a,e < M, 0<i<M. (66) 

Even if not explicitly stated, all the constants appearing in the estimates ((67)) , ([68)1 , (j69p , ([72]) , ([73]) and ((77)) 
of the following lemmas will depend also on m and Af . 

Lemma 24. For any T > and c > 1 there exists a positive constant C such tha^ 

e^r"'«'^(i, x)<c e^^f ^^'^■^*'^'^*'^(i, x), (67) 

/or any t e {0,T], x £R and N > 0. 
Proof. For any n > we have 

TZ''''{t,x) < ^f^q,,{t,x)K''^^'Ht,x), 
V m 



®Here C^ denotes the identity operator. 
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where 

qnit,x) =exp 

A direct computation shows that 
maxq„(t,x) — exp 



2{at + nO) 



\ 



2cM(t + n) 



I s'^[a^2[e^ -\\i 



8 {cM(n A- s)~ so." nS'^) 



< exp 



T a + 2 62 - iW 



8(cM - a) 



for any t E {Q,T], n > and a, 9, £ in (|551) . Then the thesis is a straightforward consequence of the fact that 
qn{t, x) is bounded on (0, T] x R, uniformly with respect to n > and a, 9, i in ([M]) . D 

Lemma 25. For any T > 0, k E N and c > 1, there exists a positive constant C such that 



d',rr''it,x) < 



c 



{at + n9)''/^ " 



\t,x), 



for any a; G M, t g]0, T] and n > 0. 
Proof. For any fc > 1 we have 



d'.Tr^%t,x) 



{at + nO)''/^ 



\t,x)pk 



'x- (f +£ei-^)t' 
^/cxt + n9 



where pk is a polynomial of degree k. To prove the Lemma we will show that there exists a positive constant 
C, which depends only on to, M, T, c and fc, such that 



cc - (f + ^e^ - ^) i 



Vat + n^ 



rr''(t,a;)<crr='^'^(t,x), j<fc 



Proceeding as above, we set 



i^+lei -l\t 



\J at 



rr^^(i,:r)=rr^'''^(t,a;)(?„,,(t,a;), 



where 



Qn.j 1,^7 *^J 



x - (f + £ei - f) t 



\/ai + J^^* 



exp 



(a;- (f +£ef -£)tV (x - (^ + ^e"* - £) iV 



2(af + ?i6') 



2{cat + nc9) 



Then the thesis follows from the boundedness of q„.j on (0,r] x R, uniformly with respect to n > and 
a,9,£ in (|55|. Indeed the maximum of q„j can be computed explicitly and we have 



lim ( max g„ ,(i, a;) 



cj 



c-1 e 



D 
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Lemma 26. For any T > and N £ N, there exists a positive constant C such that 

UQ^V-'^'^t^x) < Cr"^2(iv+i)'''^(t,a;) 
for any t € (0, T] and x G M. 
Proof. We first prove there exists a constant Cq, which depends only on m, M, T and N , such that 



(69) 



n'\t^x)<coTr'''^'^'\t,x\ 

for any t g]0, T], x G R, n > 1 and a, 6*, £ in ^^. To prove ([701) we observe that 



(70) 
(71) 



T^^^{t,x) < 



1 



(a;-(f +^e(^+i)»-^)i)'\ , , 



where 



(x - (f + £el - l) t) (a; - (I + £e(^+i)« - ^) t) 



qnit,x) =exp 
Now it is easy to check that 

niax(7„(i, a;) = exp 



2{at + {n + N)e) 



2{at + 2n{N + 1)9) 



/fM 



2/)2 



^ ^ ' < exp ^ ^ 



V 



2{n- N + 2nN)e 



V 



2iV6' 



for any i > 0. Thus g„ is bounded on (0, T] x R, uniformly with respect to n G N and a, 9, i in (|66)) . To see 
the above bound, simply observe that 



^^^^^<yWTT). 

^/at + {N + n)9 



The proof of ([71]) is completely analogous. Finally, by ((7n)) - ((7T|) we have 

n=l 

< Co (e-''itT",'^^''+''>'^'{t,x)+ete'''f2^-^K-'^''+'^'-'\t,x) 
<Co(l + Mr)r"'2(A'+i)«'^(t,a:). 



Lemma 27. For any T > and N > 2, there exists a positive constant C such that 

e^T"^'^'{t,x) < ce2jv9r"'2^''^(i,x) 

for any t G (0, T] and x G M. 



D 



(72) 
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Proof. By ^ 



n=Q 



n=0 



Lemma 28. For any T > 0, N > I and c> 1, there exists a positive constant C such that 



N 



Vl'''\t,x)<CTT^''^\t,x) 



\C0L.cQ .i I 



u 



(73) 



^\/ at + nO , 
for any x G M, t G (0, T\ and n > 0. 

Proof. We first show that there exist three constants Ci = Ci{m, M,T,N,c), C^ — Ci{N.,c) and C3 = 
C3(m, M, T, N, c) such that 



^/at + nO 



g 2cl/3(Qr + 7i8) < (^26 2c2/3(QT + ne) 



i+^c 2 -l\t 



g 2c2/3(Qr+,ie) < (^gg 2c(Qt+ne) 

for any x G M, < G (0, T] and n > 0. In order to prove (fM|) we consider 



(74) 
(75) 

(76) 



q„{t,x) = exp - 



(x- (f +^ef -^)i) 



2{at + nO) 2ci/3(Q,i + nO) 



and show that 



maxq„(t, x) = exp —i^-ttt. r-, — r < exp 

^6R^^' ' ^ 2(ci/3 _ 1) (to + „6l) ' - ^ 



(f +£et 



2 \ 



\ 



2(ci/3 _ 1) 



/ 



for any t G (0,r]. Thus (?„ is bounded on (0,T] x M, uniformly in n > and a,9,£ in (|55]). The proof of 
([76)1 is completely analogous. Equation ([75|) comes directly by setting 

N 



C2 = max ( a^e l^m^ 2.^/^ \ = e"^ 



ael 



Vcl/3 _ 1 



Now, by (|74l) we have 



^at + n0 



Af 



r^^'''(i,x)<Ci 



TV 



2cl/3(c.T + ne) 



VmTT^J y/2Tr{aT + 



(bydZB) 



<CiC2 



g 2c2/3(QT + 7ie) 

^27r(ar + n6') 
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(bydZi) 



n 



Lemma 29. For any T > 0, c > 1 and j G N U {0} there exists a positive constant C such that 

\x\elT'^-0'\t,x) < C (e2,cer^"'^''''(i,x) + Q2U+i)ceT'''^^'''-\t,x)), (77) 

for any t G (0, T] and x G M. 
Proof. By Lemma E51 there is a constant Co, only dependent on m, M, T and c, such that 

-' (£i)' 



\QlT'^-fi-\t, X) < Coe-^* ^ SllJ_^^t + {n + j)erT^'''{t, X) 



n! 

n=0 



< CoVm(Vt + j)e^^gr="'^^'^(i,x) + CoM*Te^Jir^"'^«'^(t,x), 

for any t G (0,T] and x G M and a,9,£ in (pS|) . Therefore, the thesis fohows from Lemma [^ for j = and 
from Lemma [^7] for J > 1 • D 

Lemma 30. For any T > and N,k > I we have 

e^f"'^^\t,x) < VkTTe^+''t''''^'\t,x), t G]o,r], x g m. 

Proof A direct computation shows that 



r:+jv(^,:^) Vat + {n + N + k)e ^ r-— 

max - — a — , — < V K + 1, 

-e« KiN+k (i> ^) ^at + {n + N)e 

for any t<r, n>0, Af>l and a, 6*, £ in (|66|) . This concludes the proof. D 

Proposition 31. The solution of ((27|) with n — 1, Ai as in (fT4 |) - ([T3| anrf h — Sy, takes the form 

vi(t,x-T,y)= (at - t){x - y) + ^^ ^ ^^ j\ Aip^°\t,x;T,y), 

where J is the operator 

J = ao{2d^ - 1) - Ao (e^ - 1 + S^d^es2^ , (78) 

and 6^2 is the convolution operator defined in (|65p. 

Proof. Under the assumptions in this appendix, (/>„ in (J13p takes the form 

MO = an(-C' - *0 - «^A„ (e^ - l) + A„ (er^ - l) , 



(79) 
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with 

n! n\ 

Now, by (|29p with n = 1, we have 



vi (t, C; T,y) = J e(^-*)^°(«) {zd^ - y) (^i (e)p<") (s, C; T, 2/)) ds 



Then, recalhng that 



we get 



V 27r 



^i(t,e;r,y) = ((T- i) {id^ ~y)- '-^^^Y^^m) MOp^'\t,^;T,y). 

The Proposition fohows by (j79p and by inverse Fourier transforming from ^ into the original coordinate 

X. n 

Proposition 32. For any c > 1 and t > 0, f/iere exists a positive constant C, only dependent on 

c, r, rn. A/, ||Ai||oo a?^c? ||ai||oo, such that 

\{x - y)^-"{d,, - d,)vn{t,x;T,y)\ < C{1 + WXiWooGcm) f=^^'=^'=^^(t, a;;T,2/), (80) 

for any n G {0, 1}, a;, y G M and t, T G R with <T ~t<T. 

Proof. Recalhng the expression of vo{t, x; T, y) = p^^\t, x; T, y) in (|M)). the case n = directly follows from 
Lemmas [25l [28l and [24l with iV = 0. For the case ri = 1, by Proposition [31] we have 

(x - y){d^^ - a,)wi(i, x; T, y) = (T - i)(x - y)((x - y){d^. ~ d^) + 2d, - l)Aip^°\t, x- T, y) 

+ ^^^^(^ - y)J(^-- - d,)Aip^°Ht,x;T,y), 

with J as in ([75]) and ^li acting as 

yiiM(x) = ai(92;3; — dx)u{x) — Ai ( ie~ — 1 j dxu{x) — Cs2u{x) + u{x) j . 

In the computations that follow below, in order to shorten notation, we omit the dependence of t, x, T, y in 
any function. By the commutative property of the operators dx and C, and by applying Lemmas [^ [^ and 
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[5] with A^ = 1, there exists a positive constant Ci only dependent on c, r, m, Af, || Ai||oo and ||ai||oo such that 

\{T-t){x ~ v){{x ~ y)(a,, - 9,) + 29, - l)ai(9,, ~ 9,)p(o)| < Cir™^^)'^^''^^^), 
(T - i)(x - y)((a; - y){d^,^ - d,) + 2d, ^ l)Ai ((e^ - l) 9. + l) p(°)| < Ci(r-('')'=^''^(^) + T--(y)AcS^Mv)^^^ 



(81) 

for any a;, y e K and t,T ^M. with < T — t < r. Analogously, by the commutative property of dx and C, 
and by applying Lemmas [28 l [25 | [29l and [27l with N = 2, there exists a positive constant C2 only dependent 
on c, T, m, M, || Ai||oo and ||ai||oo such that 

|(T - t)ix - 2/)((x - 2/)(9,, - dx) + 2dx - i)Aie52p(o)|||Ai|U < C2 (e,,.r™(^)^^*''^('^) + e4e^2r-(^)-4-*''^(^)), 



(82) 

for any x,y eR and i, T e M with < T - 1 < r. Then, dHHl) follows from (gll) and ([Hll) by applying Lemma 
[M] with iV = and A^ = 1 respectively. D 

Proposition 33. For any c > 1 and t > 0, there exists a positive constant C , only dependent on 

c, r, m, M, ||Ai||oo a?^c? ||ai||oo, such that 

{x - yf-^' ((e* ~l)dx + es2 ~ 1) «„(i,x;r,y)| < C(l + C^m) f^'''^'^'^'^(i,x; T, y), (83) 

for any n £ {0, 1}, a;, y G M and t, T G M mi/i < T - t < t. 

Proof. For simplicity we only prove the thesis for n = 0. The proof for n = 1 is entirely analogous to that 
of Proposition 1321 Once again, hereafter we omit the dependence of t, x, T, y in any function we consider. 
Recalling the expression of va{t,x;T,y) = p^°^{t,x;T,y) in ([64]). by Lemmas [25l [28land[29l there exists a 
positive constant Ci only dependent on c, r, to, M such that 

(a; - yf ((e^ - l) 9, + 652 - l) v^ < Ci(t^'^(v)AcS^Mv) ^ (^ ^ eiec^Or""^^^''*^"^'^^'^'), 



for any x, y € R and i, T e M with < T - t < t. Then, ([HSD follows from Lemma [2il with iV = and with 
7V = 1. D 

Proposition 34. For any c > 1 and t > 0, i/iere exists a positive constant C , only dependent on 

c, r,m,M, ||Ai||oo,||A2||oo, Ikilloo and\\a2\\oo, such that 

C'"+HT-t) 



\Zn{t,x;T,y)\ < 



[l + \\Xi\\o.e:+/)r^^''^^''^''{t,x;T,y), 



(84) 
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for any n > 0, x, y € M and t,T £ M. with < T — t < t. 

Proof. Let us define the operators 

Lo^dt+Ao, Li^dt+Ao + {x~y)Ai. 

Let us recall that, by (pS)) and (P7| with n = 1, we have 

LqVo^O, LqVi = -{Li - Lo)vq. 

Thus, by dSS]) we have 

Zo{t,x;T,y) ^ Lp^^\t,x]T,y) ^ Lvo{t,x;T,y) + Lvi{t,x;T,y) 

= {L- Li)vo{t,x;T,y) + {L ~ LQ)vi{t,x;T,y), 

where (L ~ Lq) and (L — Li) are explicitly given by 

[L - io) = {a{x) - a{y)){d,^ - 8^) + (A(.t) - X{y)) ((e^ - l) 9. + 6^2 - l) , (85) 

(i - Li) = (a(a;) - a(j/) - a'{y){x - y)){d^x ~ d^) 

+ {\{x) - X{y) - X'{y){x - y)) ((e* - l) 9. + 6^^ - l) . 

Thus, by the Lipschitz assumptions on a, A and their first order derivatives, we obtain 

\Zo{t,x;T,y)\ < ||a2||oo|a:; - y\'^\{dxx - dx)va{t,x;T,y)\ + ||ai||oo|a; - y\\{dxx - dx)vi(t,x;T,y)\ 
+ \\X2\Ux - yp I ((e^ - 1) dx + 652 - 1) vo{t, x; T, y) 
+ \\Xi\\oo\x ~ y\ |((e^ -l)dx + 652 - 1) vi{t,x;T,y) . 

and, as ||Ai||oo — implies IIA2II00 — Oj by Propositions 1321 and 1331 there exists a positive constant C, only 
dependent on c,t, to, M,||Ai||oo,||A2||ooi l|ai|loo and ||a2||ooi such that (|84p holds for n = 0. To prove the 
general case, we proceed by induction on n. First note that, by (j26p we have 



\Lp^'Ht,x;T,y)\^\{L-Lo)p'">{t,x;T,y)\ 



(0)/ 



(and by (j85p and the Lipschitz property of a, A) 



< \\ai\Ux - y\\{dxx ~ dx)p^"Ht,^;T,y)\ 

+ \\Xi\Ux - y| I ((e^ - 1) d, + es2 - 1) p(°)(i,x; T, y) 



(and by applying Lemmas [24l 



<Co 



land [291 with N ^ 0,1) 
1 



Vt^ 



+ llAi|UecAf r 



^cM,cM,cM 



{t,x;T,y), 



(86) 
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for any .x,y G R and t,T G M with < T — t < r, and where Co is a positive constant only dependent on 
c,T,m,M, ||Ai||oo and ||ai||oo- Assume now ([M]) holds for n > 0. Then by (|63l) we obtain 

\Z„+i{t,x;T,y)\< [ / |Lp(°)(f,x; s,Oll^n(s,C; T, y)|dCds 



(and by inductive hypothesis and by (j86p ) 



<^^^ TiT-s)^ I (^^ + |lAi|Ue,M)f='^'=^'=^(i,x;s,e)- 



It JR \\/s-t 

■ (1 + ||Ai||ooe^+/) r^^''^^''^'\s,^;T,y)d^ds. 
Now, by the semigroup property 

/ egf «^«'^(i, x; s, Oe^f"'''\s, C; T, y) d^ = e,^+^f "■«-^(i, x; T,y), k,N> 0, (87) 

JR 

and by the fact thalu 

with K = V2tt, we obtain 

|Z„+i(t,.;T,,)| < ^^ f^^^7==^ (1 + l|Ai||ooe:^+/)') f^^^-^^-^^(t,.;T,,) 
Vn! \ Vn +1 J 

||Ai||oo (LcA/ + ||Ai||ooL^jVj j i ■ (t,x,l,y). 



. i\ \ n + 2 
Now, by Lemma [30] we have 



g„+if cA/,cM,cM(^^ x; T, y) < 2 e';+/f ^^^■^^^■^^'^(i, x; T, y), 

CcAf f ^^^■^^^■^*^(i, x; T, y) < V^rT2 en+2^cM,cM,cM^^^ ^. j.^ ^y 



Inserting the above results into ([55)) we obtain 

|Z„+i(t,x;r,y)|< ^ ^°^ ^^ (.t + 2HAi|U(/. + V7(l + HAi|U))e^+/)r*^'-^'^'-^'^(t,x;r,y) 



< ^"^'^;^'^ ^^f ^ (1 + ||Ai|Ue;'+/) f^*^-^-^'^(i,x;r,y), 

where 

Ci = 2Co(aj + V7(1 + |1Ai|U)). 

Now, without loss of generality we can assume C\ < C, and thus we obtain (|84| for n + 1. D 



^Here Te represents the Euler Gamma function. 
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We are now in position to prove Theorem [Til for A^ = 1. Indeed, by equations (pO]) . (|FT|) and Proposition 
]we have 

\pit,x;T,y)-p<-'Ht,x;T,y)\ 
„_n Vn\ Jt Jr 



n=0 



(and by Lemma [24] with iV = and A^ = 1 respectively) 

„=o v«' Jt Jk 

(and by §7^) 

= 2{T-t) [e ^"^'^~ '^" (1 + HAilUe^+Z) r-^^--^^-^'^(t,x;r,^)') , 

for any x, y S R and t, T G M with < T - 1 < t. Since 

n— 

can be easily checked to be convergent, this concludes the proof of Theorem [Til for iV = 1. The proof for 
the general case is based on the same arguments. However, in the general case the technical details become 
significantly more complicated. Therefore, for sake of simplicity, we present here only a sketch with the 
main steps of the proof. We repeat the same iterative construction of p{t, x; T, y), but using our A^-th order 
approximation as a starting point. Thus, we replace p^^\t,x;T,y) with p'^^\t,x;T^y) in ([TO]) , where $ is 
now defined as 

oo 

'^{t,x-T,y) ^Y.^n(t,x-T,y), 

n=0 

where 

Z^{t,x;T,y):^Lp^^'>{t,x;T,y), 

Z^_^,{t,x;T,y):^ f f Lp^^\t,x;s,Oz!!{s,e,T,y)d^ds. 

Jt Jr 

Now, proceeding by induction, one can extend Propositions [32] and Proposition [33] to a general n G N. 
Eventually, after proving the identity 

N 

Lp^^^ {t, X- T,y) = Y,{L~ Lr,)v^«-^^ (t, x; T, y) + {L - Lo)v^'^ (t, x; T, y), 

one will be able to prove the estimate ([5^ for \Z^{t,x;T,y)\, from which Theorem [HI would follow exactly 
as in the case A^ = 1. 

D 

36 



D Proof of Equations ([42]) and (l43]) 



Below, we will use the following repeatedly 

0„(^)e*«- = a,M - 9,)e*«^ 0'„(Oe'^" = ^M'^d. - l)e^«^ 



K(Oe'^" = -2a„e'«^ 



Note that the above relations hold only when (j)n{£,) = an(— ^^ ^ *0j ^^ i^ the case when X is a diffusion 
with not jumps and no possibility of default. Now, using the above relations, we find 

1 



ui{t,x) 



2tt 

1 

2^ 






/2tt 
1 



27r jR Jo 
^ / ddt{x -x)- i*tVf,(C))0i(Oe^^"+**««^/i(e) 

V^TT JR ^ ^ 

(^i(a;-x) + ^t'^ao{2dx - l)jai(5^ - 5a;)Mo(i,a;), 
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which estabhshes (|42l) . Likewise, for U2 we have 

V ^TT Jv 



ZTT JR JO 

"t JR Jo 

^ f dC f ds f dre'^''+^'-''^'''°^^\id^ - x)MOe^"'''^'''"^^\id^ - x)MOe'''''°^^'^h{0 
2tt jr Jo jo 

]= [ d£, [ dsMOe"^°^^'>HOi-idi ~ xfe'^'=+^'-'^'^°'^^^ 

"t jR Jo 

27r jR Jo JO 
\J l-K Jr \ 

- \^t^ {x - x) 01(0 - ii^0;,(O0'i(e) - \t^Mi)'t>m)MOe'^"'^"^°''^m 

\2 I j-2/^ -^„ /o« 1^ I l+3„2/osi i\2 , i2„ \ left 



^\t{x- xY + t''(a; - x)aQ{2d^ - 1) + ^^■^0^(29:, - 1)^ + t^a^ 1 a2(9^ - a^)uo(i, x) 

+ ("ii^Cx - xfa^{dl - d.,) + it^ {x ~ x) aiidl - a,)ao(2a, - 1) + ^t^aiid^ ~ d,)al{2d., - 1)^ 
+ k^ {x-x)ai{2d^ - 1) + ii3ao(2a, - l)ai(25, - 1) + k'ai(9^ - a,)ao )ai(5' - d^)uo{t,x), 



which estabhshes (|43l) . 
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Figure 1: A comparison of the exact (solid), 2nd order two-point Taylor series approximation (dashed), 4th 
order Taylor series approximation (dotted), and 4th order Hermite polynomial expansion of the function 
exp(— x). For the two-point Taylor series approximation, we expand about xi = — 1 and X2 = 1. For the 
(usual) Taylor series and the Hermite polynomial approximations we expand about x ~ 0. 
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Figure 2: Using the model considered in Section [7.1.1l we plot p^") (t, x; T, y) (solid black) and p("^^) (i, x; T, y) 
(dashed black) as a function of y for n = {1, 2, 3, 4} and t = {1.0, 3.0, 5.0} years. For all plots we use the Taylor 
series expansion of ExampleEl Note that as n increases p*^"-* and p("~^) become nearly indistinguishable. In 
these plots we set /? = 0.7 All other parameter values are those listed in Section [7. 1.1 1 
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Figure 3: Implied volatility (IV) is plotted as a function of log-strike k := \ogK for the CEV-like model with 
Gaussian- type jumps of Section [7.1.11 The solid lines corresponds to the IV induced by u'-'^'>{t,x), which is 
computed using the one-point Taylor expansion (see Example [6]) . The dashed lines corresponds to the IV 
induced by u^'~''{t,x) (again, using the usual one-point Taylor series expansion). The crosses correspond to 
the IV induced by u'-*^'^-'(i,x), which is the price obtained from the Monte Carlo simulation. 



43 



t 


k 


^(3) 


u MC-95% c.i. 


IV[u(3)] 


IV MC-95% c.i. 




-0.6931 


0.0006 


0.0006 - 


0.0007 


0.5864 


0.5856 


- 0.5901 




-0.4185 


0.0024 


0.0024 - 


0.0025 


0.4563 


0.4553 


- 0.4583 


0.2500 


-0.1438 


0.0111 


0.0110- 


0.0112 


0.2875 


0.2865 


- 0.2883 




0.1308 


0.1511 


0.1508- 


0.1513 


0.2595 


0.2573 


- 0.2608 




0.4055 


0.5028 


0.5024 - 


0.5030 


0.4238 


0.4152 


- 0.4288 




-1.2040 


0.0009 


0.0009 - 


0.0010 


0.5115 


0.5176 


- 0.5210 




-0.7297 


0.0046 


0.0047 - 


0.0048 


0.4174 


0.4178 


- 0.4199 


1.0000 


-0.2554 


0.0314 


0.0313- 


0.0316 


0.3109 


0.3102 


-0.3117 




0.2189 


0.2781 


0.2775 - 


0.2784 


0.2638 


0.2620 


- 0.2649 




0.6931 


1.0034 


1.0030 - 


1.0041 


0.3358 


0.3296 


- 0.3459 




-1.3863 


0.0074 


0.0081 - 


0.0083 


0.4758 


0.4851 


- 0.4870 




-0.8664 


0.0224 


0.0224 - 


0.0227 


0.4031 


0.4029 


- 0.4045 


3.0000 


-0.3466 


0.0776 


0.0773 - 


0.0779 


0.3280 


0.3274 


- 0.3288 




0.1733 


0.3097 


0.3094 - 


0.3107 


0.2690 


0.2685 


- 0.2703 




0.6931 


1.0155 


1.0150- 


1.0169 


0.2558 


0.2540 


- 0.2604 




-1.6094 


0.0160 


0.0164 - 


0.0166 


0.5082 


0.5111 


- 0.5128 




-0.9324 


0.0439 


0.0436 - 


0.0440 


0.4118 


0.4107 


- 0.4121 


5.0000 


-0.2554 


0.1504 


0.1497- 


0.1507 


0.3203 


0.3194 


- 0.3208 




0.4216 


0.6139 


0.6123- 


0.6142 


0.2521 


0.2500 


- 0.2524 




1.0986 


2.0050 


2.0032 - 


2.0057 


0.2297 


0.2163 


- 0.2342 



Table 1: Prices (u) and Implied volatility (IV[m]) as a function of time to maturity t and log-strike k := log K 
for the CEV-like model with Gaussian-type jumps of Section [7. 1.11 The approximate price w^^^ is computed 
using the (usual) one-point Taylor expansion (see Example[5l For comparison, we provide the 95% confidence 
intervals for prices and implied volatilities, which we obtain from the Monte Carlo simulation. 
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Figure 4: Implied volatility (IV) is plotted as a function of log-strike k :— \ogK for the CEV-like model with 
Variance Gamma-type jumps of Section [7.1.21 The solid lines corresponds to the IV induced by u^^^(t, x), 
which is computed using the two-point Taylor expansion (see Example [6]) . The dashed lines corresponds 
to the IV induced by u'-^>{t,x) (again, computed using the two-point Taylor series expansion). The crosses 
correspond to the IV induced by u'^^'^\t,x), which is the price obtained from the Monte Carlo simulation. 



t 


k 


U(2) 


u MC-95% c.i. 


IV[m(2)] 


IV MC-95% c.i. 




-0.6931 


0.0014 


0.0014 - 0.0015 


0.4631 


0.4624 - 0.4652 




-0.4185 


0.0070 


0.0070 - 0.0071 


0.4000 


0.3995 - 0.4014 


0.5000 


-0.1438 


0.0363 


0.0362 - 0.0365 


0.3336 


0.3331 - 0.3346 




0.1308 


0.1702 


0.1697- 0.1704 


0.2727 


0.2707 - 0.2736 




0.4055 


0.5011 


0.5004 - 0.5012 


0.2615 


0.2291 - 0.2646 




-0.9163 


0.0028 


0.0027 - 0.0028 


0.4687 


0.4678 - 0.4702 




-0.5697 


0.0109 


0.0109- 0.0110 


0.4057 


0.4050 - 0.4068 


1.0000 


-0.2231 


0.0473 


0.0472 - 0.0476 


0.3434 


0.3428 - 0.3444 




0.1234 


0.1970 


0.1965 - 0.1974 


0.2836 


0.2825 - 0.2847 




0.4700 


0.6033 


0.6025 - 0.6037 


0.2452 


0.2355 - 0.2506 



Table 2: Prices (m), Implied volatilities (IV[w]) and the corresponding confidence intervals from Figure H] 
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Figure 5: Implied volatility (IV) is plotted as a function of log-strike k :— logiiT for the model of Section [7?2] 
The dashed line corresponds to the IV induced by u'^^'>(t, x). The solid line corresponds to the IV induced by 
u'^^\t,x). To compute u'^^\t,x), i € {0,2}, we use the two-point Taylor series expansion of Example [6] The 
crosses correspond to the IV induced by the exact price, which is computed my truncating (1501) at n — 8. 
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Figure 6: Left: survival probabilities u{T — t,x) :— Qx[C > T\(^ > t] for the JDCEV model described 
in Section mi The dotted line, dashed line and solid line correspond to the approximations w^°)(T — t,x), 
u^^\T—t,x) and u'^'(T— i,a;) respectively, all of which are computed using CoroUarvfTOl The crosses indicate 
the exact survival probability, computed by truncating equation (|54p at n = 70. Right: the corresponding 
yields y(")(i, x; T) :— — log(w("^(T — i, x))/{T — t) on a defaultable bond. The parameters used in the plot 
are as follows: x — log(l), /3 = —1/3, h = 0.01, c = 2 and a = 0.3. 



T-t 


Y 


Y - r(o) 


r-y(i) 


y_y(2) 


1.0 


0.1835 


-0.0065 


0.0022 


0.0001 


2.0 


0.1777 


-0.0123 


0.0048 


0.0003 


3.0 


0.1720 


-0.0180 


0.0071 


0.0003 


4.0 


0.1663 


-0.0237 


0.0089 


-0.0001 


5.0 


0.1605 


-0.0295 


0.0099 


-0.0006 


6.0 


0.1548 


-0.0352 


0.0102 


-0.0011 


7.0 


0.1493 


-0.0407 


0.0101 


-0.0013 


8.0 


0.1442 


-0.0458 


0.0095 


-0.0011 


9.0 


0.1394 


-0.0506 


0.0087 


-0.0005 


10.0 


0.1351 


-0.0549 


0.0077 


0.0007 



Table 3: The yields Y{t,x;T) on the defaultable bond described in Section TTM exact {¥) and nth order 
approximation (F^")). We use the following parameters: x — log(l), (3 — —1/3, b — 0.01, c = 2 and 6 = 0.3. 
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Figure 7: Top: imphed volatihty for the CEV model for maturities t = 0.5 and t = 1.0. The solid line is the 
exact implied volatility a, computed by truncating Q in equation ()55p at n = 100 terms and then inverting 
Black-Scholes numerically. The dashed line is the 2nd order implied volatility expansion u"' , computed using 
equation (|44p . The dotted line is the Hagan- Woodward expansion of implied volatility a^^ , computed using 
equation ([56]) . Bottom: relative errors of the 2nd order (dashed) and Hagan- Woodward (dotted) implied 
volatility expansions. In all plots we use the following parameters: /? — 0.1, S — 0.3, Xq — log(l.O). Units of 
the horizontal axis are log-strike: logiiT. 
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Figure 8: Top: imphed volatihty for the CEV model for maturities t = 2.0 and t = 3.0. The solid line is the 
exact implied volatility a, computed by truncating Q in equation ()55p at n = 100 terms and then inverting 
Black-Scholes numerically. The dashed line is the 2nd order implied volatility expansion u"' , computed using 
equation (|44p . The dotted line is the Hagan- Woodward expansion of implied volatility a^^ , computed using 
equation ([56]) . Bottom: relative errors of the 2nd order (dashed) and Hagan- Woodward (dotted) implied 
volatility expansions. In all plots we use the following parameters: /? — 0.1, S — 0.3, Xq — log(l.O). Units of 
the horizontal axis are log-strike: logiiT. 
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